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(54) Method and system of compensation of thickness of an organ 



(57) Method and system for compensating for the 
thickness of an organ in a radiology instrument, in which 
an image of the radiological thicknesses of the organ 
through which the X-ray beam has passed is calculated 
(15) on the basis of a digitized image, the thickness im- 
age is filtered using a low-pass filter (1 6) in order to ob- 
tain a low-frequency image, the low-frequency image is 
subtracted (1 8) from the radiological thickness image in 
order to obtain a contrast image, the low-frequency im- 



age is processed (1 7) using a pre-recorded table taking 
into account a contract % selected by a user in order to 
obtain an image with reduced dynamic range, and the 
image with reduced dynamic range is added (1 9) to the 
contrast image in order to obtain a compensated thick- 
ness image, the pixels having a level below or above a 
predetermined threshold being returned at least to the 
value of the said threshold, while preserving the differ- 
ences and real ratios between the anatomical struc- 
tures. 
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Description 

[0001 ] The present invention concerns the field of ra- 
diological imaging to visualize an organ or part of an or- 
gan, generally of the human body. 
[0002] Radiography is conventionally carried out with 
sensitive films exposed by X-rays after they cross the 
organs to be studied. Radiologists have been trained in 
the interpretation of such Images. The new imaging 
technologies - solid state detector and digital acquisition 
system - must be adapted to common practices and 
must provide an equivalent perception of the pertinent 
information that radiologists are accustomed to exam- 
ining. In particular, one of the requirements that digital 
systems must satisfy consists of reducing the extent of 
the gray level dynamics in order to simulate a conven- 
tional film as faithfully as possible. For this purpose, the 
digital image is displayed on a screen that the radiologist 
adjusts interactively in order to identify all the clinical 
signs by perceiving the relations between the different 
image components. But the precise perception of den- 
sity information by means of the image displayed on the 
screen is limited by the dynamics of the screen. Now, 
one must pass automatically from the high contrast of 
the original Image, for example, in the order of 30 to 50, 
to the low contrast offered by a video screen. 
[0003] Physicians commonly apply configuration 
techniques to patients with a view to limiting the dynam- 
ics of acquired images, for example, by compression of 
thick regions and/or by addition of absorbent substanc- 
es in order to compensate for thin zones. In the case of 
X-ray mammography, the breast is compressed to the 
smallest and most constant possible thickness. In the 
field of cardiology, contour filters are used to avoid prob- 
lems associated with weakly absorbent zones of the 
chest, such as the lungs. These fitters consist of plates 
of shape complementing that of the heart, made of ma- 
terials of given X-ray absorption coefficients. 
[0004] However, these techniques prove insufficient 
for digital image acquisition and processing and cum- 
bersome to use, if not uncomfortable for the patient. 
[0005] The present invention is intended to remedy 
the above-mentioned problems. 
[0006] The embodiments of the present invention 
compensate for variations of thickness, in particular at 
the limit between high-density regions and low-density 
regions of the organ X-rayed, in an adaptive manner, 
depending on characteristics of the image visualized, 
chosen by a user. 

[0007] The present invention is intended to offer an 
image on which tissues presenting absorptions different 
from one another are naturally and exploitably repre- 
sented. 

[0008] The method of compensation of thickness of 
an organ is designed for an X-ray apparatus of the type 
comprising an X-ray source and a means of detection 
of the X-ray beam after it has crossed the organ, the 
means of detection being capable of converting the X- 



ray beam into a digital electronic signal. From a digitized 
image, an image of the radiological thicknesses of or- 
gans crossed by the X-ray beam is calculated, the image 
of radiological thicknesses is filtered by a low-pass filter 

5 in order to obtain a low-frequency image, the low-fre- 
quency image is subtracted from the image of the radi- 
ological thicknesses in order to obtain a contrast image, 
the low-frequency image is processed by a prerecorded 
table taking into account a contrast % chosen by a user 

10 in order to obtain an image with reduced dynamics, and 
the image with reduced dynamics and the contrast im- 
age are added to obtain a compensated thickness im- 
age. The pixels of level lower or higher than a predeter- 
mined threshold are brought back at least to the value 

15 of the threshold, while preserving the differences and 
the real relations between the anatomical structures. 
[0009] The thickness of an organ measured by X-rays 
is called radiological thickness, in other words, taking 
into account the absorption of the materials crossed. For 

20 example, 1 cm of bone has the same radiological thick- 
ness as 1 0 cm of water. 

[0010] The thickness image can be obtained by 
means of Lambert's law: I = ^e-v* with I the number of 
photons received at a given point of the means of de- 

25 tection, l Q the number of photons which would be re- 
ceived at a given point of the means of detection, if the 
organ was not present in the field of view, u. the coeffi- 
cient of linear attenuation of the X-rays by the material 
crossed, and t the thickness of material crossed, from 

30 which one deduces: 

[0011] In l 0 - 1 n I = ui, the product u,t corresponding 
to the gray level for a pixel of a radiological thickness 
image. 

[0012] Calculation of the image with reduced dynam- 

35 ics advantageously takes into account a range of gray 
levels of width and center indicated by a user. 
[0013] In an embodiment of the invention, the low-fre- 
quency image processing for obtaining the image with 
reduced dynamics is carried out as a function of the 

40 width and center of the passband. 

[0014] In an embodiment of the invention, the low-fre- 
quency image processing for obtaining the image with 
reduced dynamics is carried out by means of a digital 
table or an analytical law. 

45 [0015] In an embodiment of the invention, the low-fre- 
quency image processing follows a monotone law. 
[0016] In an embodiment of the invention, the low-fre- 
quency image processing follows a linear law of slope a. 
[001 7] The slope a advantageously evolves in a man- 

50 ner inversely proportional to the contrast %. 

[0018] In an embodiment of the invention, the low-fre- 
quency image is stored in a memory and, on a change 
of contrast, the low-frequency image is read in the mem- 
ory and the necessary processings and calculations are 

55 carried out. 

[0019] In another embodiment of the invention, on a 
change of contrast, the image of radiological thickness- 
es is filtered by a low-pass filter in order to obtain a low- 
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frequency image and the necessary processings and 
calculations ar carried out. 

[0020] Th present invention also concerns a system 
of compensation of thickness of an organ in an X-ray 
apparatus. The X-ray apparatus is of the type compris- 
ing an X-ray source and a means of detection of the X- 
ray beam after it has crossed the organ, the means of 
detection being capable of converting the X-ray beam 
into a digital electronic signal. 
The system includes a means of calculating, from a dig- 
itized image, an image of the radiological thicknesses 
of organs crossed by the X-ray beam, a low-pass filter 
for obtaining radiological thicknesses from the image, a 
low-frequency image, a means for subtracting the low- 
frequency image from the radiological thicknesses in or- 
der to obtain a contrast image, a means for processing 
the low-frequency image according to a prerecorded ta- 
ble taking into account a contrast % chosen by a user in 
order to obtain an image with reduced dynamics and a 
means for adding the image with reduced dynamics and 
the contrast image in order to obtain a compensated 
thickness image, the pixels of level lower or higher than 
a predetermined threshold having been brought back at 
least to the value of the threshold, while preserving the 
differences and the real relations between the anatom- 
ical structures. 

[0021] In an embodiment of the invention, the low- 
pass filter is calculated to eliminate the pixels corre- 
sponding to an organ loaded with contrast medium in 
the compensation image and thus maintain them in the 
compensated image. 

[0022] Thus, the invention offers an image processing 
supplying an image equivalent to that which would be 
obtained by using an absorbent liquid on the edges of 
the organ over a part of their height and makes it pos- 
sible to simulate a physical phenomenon, which pro- 
vides a better understanding of the optimal adjustment 
of the parameters that has to be made. This image 
processing confers a natural and pleasing appearance 
to the different tissues of the organ studied. Each of 
those tissues is seen at the same time and on the same 
image with, on the one hand, a natural appearance and, 
on the other, a precision and quality enabling the user 
to derive important information from it, as if the user 
were interested only in one particular tissue. In the field 
of mammography, the glandular zones and the adipose 
zones are observed simultaneously with the possibility 
of deriving information from the same image on both of 
the zones. 

[0023] The method is adaptable to existing X-ray ma- 
chines and can be applied to the radiography of any or- 
gan at all. 

[0024] Whateverthe extent and centering of the range 
of gray levels chosen by the user for visualization, a 
choice made to favor certain types of organs or tissues 
in the organ studied, the image is processed in a mann r 
adapted to that choice. 

[0025] The present invention will be better understood 



by studying th detailed description of an embodiment 
taken by way of non-limitative example and illustrated 
by the attached drawings, in which: 

5 Figure 1 is a schematic representation of an organ 
composed of zones of different thicknesses; 

Figure 2 is a schematic representation of an organ 
composed of zones of different compensated thick- 
10 nesses; 

Figure 3 is a schematic view of a breast com- 
pressed for a mammogram; and 

15 Figure 4 is a schematic view of a system according 
to an embodiment of the invention. 

[0026] As can be seen in Figure 1 , the dynamics of 
an image of an organ can exceed the dynamics of thick- 

20 nesses representable by the display screen and repre- 
sented by the two horizontal dotted lines. Thus, blocks 
1 to 3 are at too high a gray level, while blocks 8 and 9 
are at too low a level. The physical phenomena or parts 
of the organ found in such zones are therefore not sat- 

25 isfactority displayed. By maintaining the same screen 
dynamics, it is conceivable, first of all, to display blocks 
1 to 4, then, secondly and after modification of the width 
and center of the range of gray levels of the screen, 
blocks 5 to 9. By taking the example of the mammogram, 

30 a first image would be obtained, on which only the con- 
tour of the breast would be seen, without viewing its in- 
ternal tissues. The second image would show solely its 
internal tissues, adipose zone or glandular zone, without 
its contour being perceived, which is hardly practical and 

35 complicates the work of the radiologist. 

[0027] In an embodiment of the invention, the digital 
compensation of the crude image is based on the fol- 
lowing principle: the image compensation process sim- 
ulates the addition or removal of an appropriate quantity 

40 of material in boundary regions, for example, between 
zones of high absorption and zones of low absorption 
of X-rays, so that the dynamics are reduced, while taking 
into account the constraint of preservation of the differ- 
ence and real relations between the anatomical struc- 

45 tures. Thus, it can be seen that in Figure 2, blocks 1 to 
9 are contained within the limits of the screen dynamics. 
Blocks 1 to 3, instead of being white-saturated, are 
brought back just below the upper limit of saturation and 
preserve their original relation, namely: block 1 lighter 

50 than the block 2 lighter than the block 3, etc. The same 
is true for blocks 8 and 9 which, instead of being black- 
saturated, are situated in proximity to the lower limit of 
the dynamics, so as to be visible and, while preserving 
their original difference, block 7 lighter than the block 8 

55 lighter than the block 9. 

[0028] In Figur 3, a breast 10 can be seen, com- 
pressed between a plate 11 and a table 12 which form 
part of an X-ray apparatus, not represented. The breast 
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is subjected to an X-ray beam 13 emitted by a source 
not represented. A digital detector, not represented, is 
placed on the path of the X-ray beam 13 after crossing 
the breast 10. It can be observed that some X-rays, ref- 
erenced 14, cross only a part of the thickness of the 
breast owing to its rounded shape at the tip. Conse- 
quently, these X-rays are less attenuated than those 
having crossed the total thickness of the breast, which 
risks leading to a saturation of the image portion corre- 
sponding to rays 1 4 and to the display of a black screen 
portion. In fact, a white portion corresponds to a thick 
zone and a black portion corresponds to a thin zone by 
convention. 

[0029] As can be seen in Figure 4, an application of 
the system according to the invention consists of a 
mathematical processing module 15 capable of apply- 
ing to an image received in input from a digital detector, 
not represented, a logarithm function. In fact, the level 
of each pixel of the input image is representative of the 
intensity of the X-radiation received. The logarithm func- 
tion makes it possible to pass from radiological intensi- 
ties to thicknesses. 

[0030] For a given point of the image, the number I of 
photons received by a given point of the means of de- 
tection is determined by the following equation: I = l max 
x e - ^, t being the thickness of the organ crossed by the 
X-rays and u. being the density of the organ crossed. 
Applying the logarithm function leads to u. x t = In (l max ) 
- In (I), l max being known and corresponding to the 
number of photons which would be received on a given 
point of the means of detection if the organ was not 
present in the field of view. The product u. x t is therefore 
known, being called densitometric or radiological thick- 
ness. 

[0031 ] I n practice, it may be preferred to apply a slight- 
ly different equation, based on quantity G, which is the 
gray level of a given point of the means of detection: (it 
= K x 1n (G max +.1) = K x 1n (G + 1). Gain K is applied 
to protect the variation of gray level higher than 1% of 
maximum intensity, such as K x 1n (G max + 1) - K x 1n 
(0.99 G max + 1)^1. 

The application of that logarithm function makes it pos- 
sible to obtain a histogram of thickness. 
[0032] The processing system further includes a filter 
1 6 of low-pass type for receiving the pixels X of an input 
thickness image and supplying on output pixels filtered 
or of input mask M x , an operator 1 7 placed downstream 
from the filter 1 6 and performing a function noted O, and 
a summator-subtractor 1 8 receiving the pixels of the in- 
put thickness image on the output of the operator 1 5 and 
also receiving the pixels filtered on the output of the filter 
16 in order to subtract the filtered pixels from the pixels 
of the input thickness image and to supply a contrast 
image of pixels C = X - M x . The filter 16 may preserve 
the structures of size greater than a nominal value, be- 
ing associable with th size of structures of interest 
present in the organ, for example, equal to 2 cm. 
[0033] The processing system includes a summator- 



adder 19 placed downstream from the operator 17 and 
from the summator-subtractor 1 8 receiving the contrast 
pixels C and the pixels M y . The output pixels of the op- 
erator 1 7 are not d M y with M y = <D(M X ). The summator- 
5 adder 1 9 performs the operation of addition of the con- 
trast pixels C and output pixels M y of the operator 17 
and supplies on output the output pixels noted Y, with Y 
= C + M y or even Y = X - M x + <D(M X ). 
[0034] An image transformation is thus made. The 
10 function O of the operator 1 7 is a monotone function with 
adjustable slope in order to simulate the addition or re- 
moval of an appropriate quantity of material in different 
regions, so that the dynamics of the image will be re- 
duced, while preserving the real relations between the 
15 anatomical structures. 

The preservation of the relations is associated with the 
monotony of the function O. Preferably, <D is a linear 
function of slope a, the slope a being adapted to the 
contrast chosen by the user. If a = 1 , then the dynamics 
are preserved and M y = M x . The operator 1 7 can contain 
a prerecorded table which takes into account the con- 
trast x of the image having to be displayed, as chosen 
by a user. As a variant, the operator 17 can follow an 
analytical law. 

[0035] The contrast % of the image having to be dis- 
played may be determined by a user by choosing a 
range of gray levels W w of the image and a center of the 
range of gray levels, for example, by keyboard com- 
mand or by means of a mouse, not represented. The 
slope a established by a table or by an analytical law 
evolves in a manner inversely proportional to the con- 
trast x and, more precisely, to the width W w of the range 
of gray levels. In other words, when the user changes 
the width W w of the range of gray levels from W w to W' w , 
the slope evolves from a to a', with W w /W w = a'/a. The 
general appearance of the image is thus preserved and 
so is that of the tissues observed, whatever their type, 
for example, pectoral, glandular, adipose, subcutane- 
ous or cutaneous in mammography. 
[0036] On such change of contrast, the different stag- 
es can be repeated from the thickness image stored in 
a memory, not represented, of the system which is di- 
rectly accessed or by means of a data bus. 
[0037] The filtered image can also be stored in a 
memory and repeat only the stages subsequent to fil- 
tering, which reduces the times of calculation, but slight- 
ly complicates the management of memory access. 
[0038] If a particular pixel X n is of gray level represent- 
ative of a thickness greater than a pixel X m and less than 
a pixei X p , the same relation will exist between the out- 
put pixels: Y m <Y n <Y p . The dynamics of an image is thus 
reduced in a proportion which will depend on the dynam- 
ics of an input image sent by a digital detector and the 
dynamics of an output image, so that the display means 
are capable of making it available to a user. 
The different processings are carried out on so-called 
"thickness" images, that is, in which the value of each 
pixel is representative of the thickness of the tissues 
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crossed by the X-rays. 

[0039] The output pixels Y ar sent to a display tab! 
20 from to which a user can choose th contrast % of th 
image visualized. 

[0040] The processing method is available for visual- 
izing on the same image tissues of very different radio- 
logical characteristics, by breaking free, at least In part, 
of the constraints of external luminosity associated with 
illumination of the site on which visualization takes 
place. 

[0041] The image obtained according to disclosed 
embodiments presents no artifact and is of normal ap- 
pearance. The method can be used in various fields of 
radiology, while being adapted to digital detection and 
processing, without changing the users' practices. An 
increase of contrast is prevented from being manifested 
by a reduction of dynamics of the image displayed, as 
commonly occurs. An increase of contrast is also pre- 
vented from being manifested by a saturation of certain 
zones of the image representing tissues that it is desir- 
able to observe at the same time as other tissues rep- 
resented by nonsaturated zones, as, for example, in 
mammography, the saturation of adipose zones becom- 
ing too dark if the contrast is increased on the glandular 
zones or, conversely, the saturation of the adipose 
zones becoming too dark if the contrast is increased on 
the glandular zones or, conversely, saturation of the 
glandular zones becoming too light, if the contrast on 
the adipose zones is increased. This is due to the fact 
that contrast is taken into account in image processing 
for display. 



Claims 

1 . A method of compensation of thickness of an object 
in a radiological imaging apparatus having radiation 
source and a means of detection of the radiation 
source after it has crossed the object, the means of 
detection being capable of convertfng the radiation 
source into a digital electronic signal, in which, from 
a digitized image, an image of the radiological thick- 
nesses of the object crossed by the radiation source 
is calculated, the image of radiological thicknesses 
is filtered by a low-pass filter in order to obtain a low- 
frequency image, the low-frequency image of radi- 
ological thicknesses is subtracted in order to obtain 
a contrast image, the low-frequency image is proc- 
essed by a prerecorded table or calculated in real 
time taking into account a contrast % chosen by a 
user in order to obtain an image with reduced dy- 
namics and the image with reduced dynamics and 
the contrast image are added in order to obtain a 
compensated thickness image, the pixels of level 
lower or higher than a predetermined threshold hav- 
ing been brought back at least to the value of th 
threshold, while preserving the differences and th 
real relations between the features of th object. 



2. Th method according to claim 1 , in which calcula- 
tion of the imag with reduced dynamics takes into 
account a range of gray levels of width and center 
indicated by a user. 

5 

3. The method according to claim 2, in which process- 
ing of the low-frequency image for obtaining the im- 
age with reduced dynamics is carried out In accord- 
ance with the width and center of the range of gray 

10 levels. 

4. The method according to any one of claims 1 to 3, 
in which processing of the low-frequency image for 
obtaining the image with reduced dynamics is car- 

15 ried out by means of a digital table or an analytical 
law. 

5. The method according to claim 4, in which process- 
ing of the low-frequency image follows a monotone 

20 law. 

6. . The method according to claim 5, in which process- 
ing of the low-frequency law follows a linear law of 
slope a. 

25 

7. The method according to claim 6, in which the slope 
a evolves in a manner inversely proportional to the 
contrast x- 

30 8. The method according to any one of claims 1 to 7, 
in which a low-frequency image is stored in a mem- 
ory and, on a change of contrast, the low-frequency 
image is read in the memory and the necessary 
processings and calculations are carried out. 

35 

9. The method according to any one of claims 1 to 7, 
in which, on a change of contrast, the parameters 
of the filter applied to the image of radiological thick- 
nesses by a low-pass filter are adapted in order to 

40 obtain a low-frequency image and the necessary 
processings and calculations are carried out. 

1 0. A system of compensation of thickness of an object 
in a radiological imaging apparatus, of the type 

45 comprising a radiation source and a means of de- 
tection of the radiation source after it has crossed 
the object, the means of detection being capable of 
converting the radiation source into a digital elec- 
tronic signal, including a means (1 5) for calculating, 

so from a digitized image, an image of the radiological 
thicknesses of object crossed by the radiation 
source, low-pass filter (1 6) for obtaining radiological 
thicknesses from the image, a low-frequency im- 
age, means (1 8) for subtracting the low-frequency 

55 image from th radiological thicknesses in order to 
obtain a contrast image, means (1 7) for processing 
the low-frequency image according to a prerecord- 
ed table taking into account a contrast % chosen by 
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a user in order to obtain an image with reduced dy- 
namics and means (19) for adding the image with 
reduced dynamics and the contrast image in order 
to obtain a compensated thickness image, the pix- 
els of level lower or higher than a predetermined s 
threshold having been brought back at least to the 
value of the threshold, while preserving the differ- 
ences and the real relations between the features 
of the object 



6 



EP 1 113 392 A1 




7 



EP 1 113 392 A1 



FIG.3 




8 



i 



,1 

EP 1 113 392 A1 



FIG.4 



18 I 



19 y 



My 



Mx 



20 



17 



9 



EP 1 113 392 A1 



Eur pean Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 00 31 1605 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with Indication, where 
of relevant passages 



Relevant 
to claim 



CLASSIFICATION OF THE 
APPLICATION (tntCt.7) 



EP 0 352 016 A (MINNESOTA MINING & MFG) 
24 January 1990 (1990-01-24) 

* abstract * 

* page 3, line 23 - line 36 * 

* page 5, line 31 - line 33 * 

EP 0 732 669 A (EASTMAN KODAK CO) 
18 September 1996 (1996-09-18) 

* abstract * 

* column 4, line 6 - line 54 * 

* column 5, line 8 - line 27 * 



1-10 



1-10 



G06T5/40 



Tne present search report has been drawn up for all claims 

Placaot search *~ 



TECHNICAL FIELDS 
SEARCHED (lnt,CL7) 



G06T 



THE HAGUE 



Oats o> computer of ttw aaaret- 

10 May 2001 



Gonzalez Ordonez, 0 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken alone 

Y : particularly relevant if combined with another 

document of the same category 
A . technological background 
O : non-written disclosure 
P : intermediate document 



T : theory or principle under tying the Invention 
E : eartter patent document, but published on. or 

after the filing date 
D : document cited in the application 
I : document cled lor other reasons 



A : member of the same patent tamiry, corresponding 



10 



EP1 113 392 A1 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. 



EP 00 31 1605 



The European Patent Office Is In no wayllaMe tor these particulars which are merely given tor the purpose ol Information. 

10-05-2001 



Paten I document 
cited In search report 



Publication 
date 



Patent famUy 
member(s) 



EP 0352016 



24-01-1990 



US 
CA 
DE 
DE 
IL 
JP 



4941190 A 
1309502 A 
68914206 0 
68914206 T 
90654 A 
2073477 A 



EP 0732669 



18-09-1996 



US 
JP 



5774599 A 
9018811 A 



i For more details about this annex : see Official joumaJ of the European Patent Office. No. 12/82 



Publication 
date 



10-07-1990 
27-10-1992 

05- 05-1994 

06- 10-1994 
21-06-1992 
13-03-1990 



30-06-1998 
17-01-1997 



11 



— %1 



I. 

r 

* 

ft 



